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ABSTRACT 

We present a study of the peculiar Type la supernova 2001ay (SN 2001ay). The 
defining features of its peculiarity are: high velocity, broad lines, and a fast rising 
light curve, combined with the slowest known rate of decline. It is one magnitude 
dimmer than would be predicted from its observed Amis, and shows broad spectral 
features. We base our analysis on detailed calculations for the explosion, light curves, 
and spectra. We demonstrate that consistency is key for both validating the models 
and probing the underlying physics. We show that this SN can be understood within 
the physics underlying the Amis relation, and in the framework of pulsating delayed 
detonation models originating from a Chandrasekhar mass, M(jy^, white dwarf, but with 
a progenitor core composed of 80% carbon. We suggest a possible scenario for stellar 
evolution which leads to such a progenitor. We show that the unusual light curve decline 
can be understood with the same physics as has been used to understand the Amis- 
relation for normal SNe la. The decline relation can be explained by a combination of 
the temperature dependence of the opacity and excess or deficit of the peak luminosity, 
a, measured relative to the instantaneous rate of radiative decay energy generation. 
What differentiates SN 2001ay from normal SNe la is a higher explosion energy which 
leads to a shift of the ^^Ni distribution towards higher velocity and a < I. This result 
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is responsible for the fast rise and slow decline. We define a class of SN 2001ay-like 
SNe la, which will show an anti-Phillips relation. 

Subject headings: Supernovae: individual: SN 2001ay 



1. Introduction 



While some progress has been made in the understanding of the Type la supernova (SN la) phe- 
nomenon in recent years, there does not y et exist an agreed upon s tandard model of the supernova 
explosion that can explain norma l SNe la feranch et al. 2005 . 2009 ). For cosmology, the br ightness 
decline relation plays a key role (jPhillipsI Il993l : iPhillips et al.l Il999l : iGoldhaber et al.l l200ll ) . Prom 
theory, AmLt^ is well und erstood: the light curves (LCs) are powered by radioactive decay of ^^Ni 
(jColgate fc McKeJl969l ;i. More ^'^Ni increases the luminosity and causes the envelopes to be hotter. 
Higher temperature means higher opacity and , thus, longer diffusion time scales and slower decline 



rates after maximum light (iHoflich et al.lll996l : iNugent et al.lll997l : lUmeda et al.lll999l : iKasen et al 



20091 ). The existence of a Amis-relation hol ds for virtually all scenarios as long as there is an 
excess amount of stored energy to be released (jHofiich et al.lll996l ). Although, the tightness of the 
relation can be understood within the framework of the single de generate scenario and spherical 



models (Hoflich et al. 



1996 



Hoflich et al 



20021 : 



Hoflich et al 



2010l). it falls ap art when taking into 



account burning instabilities during the deflagration phase (jKasen et al.ll2009l ). This difficulty and 
the observation of a set of extremely bright SNe la may ler id support for double degenerate sce- 
narios with progenitors well above the Chandrasekhar mass (jHowell et al.ll2006l : IScalzo et al.ll20ld : 
Taubenberger et al.ll201ll : iHowelll l201ll ). We note, however, that the inferred brightness depends 
on a unique relation between the ^^Ni mass Mjqi, and the intrinsic color: B — V at maximum light. 
At least in a few cases, the apparent b rightness can be uri derstood within the framework of Mch 
mass WDs with intrinsically red color (jQuimbv et al.ll2006l ). 



Additionally, some progress has been made in understanding variations among the groups, 
with sug gestions that some of the spectral diversity is due to metallicity, central dens i ty, and asym- 



metr ies (IHoflich et al.lll998l : iLentz et al.ll2000l : iHoflich et al.ll2010l : iMaeda et al.ll2010l : iMaund et al 



2010l ). The nature of the progenitor system is also controversial, with recent work on rates and 



Yuneelson & Livio 


200nl: 


MaozlboOf^ 


Ruiter et al. 


2011 




Maoz et ahlboil) 



the double degenerate scenario. (lYungelson fc Livio 



200nl: iMaozj IgOOgI : IPuiter et al.l[2009l : iMaoz et al.lboid : iMennekens et al 



1998 



201C 



Maoz et ahlboil ). Nevertheless, theoretical work continu es to favor the single 
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Despite the uncertainties in theory, most of the kno wn SNe la obey the brightness decline rela- 



tion. The light curves are self-similar within ±0.3 mag ( Goldhaber &: Perlmutter 1998j ^ 



19991 ). Deviations of this order can be expected from variations of the progenitor (jHoflich et al 



Riess et al 
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19981 : iBrachwitz et alJbood : IXhielemann et al 



works for both the local and high-z samples (jPerlmutter et al 



2003: Seitenzahl et al 



1999 



201 ll). The stretch ing method 



Riess et alJll998l ). The self- 



similarity holds even for supernovae for which super-Chandrasekhar mass progenitors have been 
suggested. 



The subject of this paper, SN 2001ay, shows that nature is even inore d iverse (IKrisciunas et al 



2011 



From the redshift of the host galaxy IC4423 (jFreedman et al.ll200ll ). the distance modulus 
IJL = 35.55 ± 0.1 is well known. The reddening in our galaxy is found to be E(B — V) = 0.026 it 
0.006 mag, and interstellar absor ption in Na I suggests low reddening in the host galaxy of about 
E(B — V) = 0.072 lb 0.008 mag (|Krisciunas et ahlboilh . Together with the observed magnitude, 
the intrinsic brightness at maximum light is My = —19.17 mag with a B — V color of —0.02 mag, 
which is comparable with typical core-normal supernovae. 

However, the light curve shape is unlike other SNe la: Its measured decl ine rate of Ator = 0.68 , 



is slower than any known SNe la, combined with a fast rise of some 16 days (jKrisciunas et al.ll2011 



Based on its Amis, SN 2001ay should be brighter than observed by about 1 mag. Moreover, the 
linewidths near maximum light put it solidly in the broad-line class of SNe la. 

We show that SN 2001ay can be understood with the same physics underlying the Amis- 
relation, a nd in the framework of parametrized pulsating delayed detonation models similar to 
SN1990N (IKhokhlov et al .11 19921 ). but with an unu sual progenitor star. We treat consistently the 
explosion, light curves, and spectra (jHofiichl 1 19951 ). We show that consistent models reproduce 
the observed light curves and spectra reasonably well. Furthermore, we show that inconsistent 
calculations lead to spectral features which would lead to rejection of the explosion model. Finally, 
we summarize our findings and discuss possible implications for the understanding of SNe la and 
cosmology. 



2. Motivation for our Model for SN 2001ay 

For our models, an understanding of t he brightness decline relation is important. The Phillips 
relation (|Phillipslll993l : iPhillips et al.lll999l ) provides an empirical link between peak brightness and 
decline of the LC after maximum. Most SNe la fall within 0.2 mag of this relation. 

As discussed in the introduction, the Amis-relation can be understood as an o pacity effect if 



energy in stored energy is available in excess of the instantaneo us radioactive decay (lH5fiich et al 



1996 



Nugent et al. 



1997 



Umeda et al.lll999l : iKasen et al.ll2009l ). The latter is a key to understand 



SN 2001ay. We will discuss that while the second condition is valid in most scenarios for SNe la, 
it is not guaranteed. 

For our discussion a useful quantity is the relation between instantaneous energy deposition 
by radioactive decay and the brightness at maximum: 



II 
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where a accounts for the fact that some of the 7-ray energy deposited prior to maximum hght can 
be stored in the thermal energy and trapped radiation energy which cannot esca pe faster than a 



diffusion time. In a typical dela yed detonation model the value of a is about 1.2 (jKhokhlov et al, 
I993I : iHofiich Khokhlovl[l996l l. 



The role of the opacity condition can be seen from Arnett's analytic solution (lArnettlll98Cl . 



19821 ). which also provides further insights. Constant opacity in a polytropic model with homoge- 
neous energy input does not produce a brightness decline relation. In this model, energy balance 
between radioactive decay and adiabatic expansion and cooling cancel exactly in a radiation domi- 
nated regime with pure geometrical diluti on. At maximum, the rat io between maximum brightness 
and instantaneous energy input is a = 1. iPinto &: EastmanI (|200d ) reconsidered Arnett's one-zone 



models and confirmed previous explanations. Declining opacit y provide extra ene rgy at maximum 
light because it accelerates the recession of the "photosphere" (jHoflich et al.lll993l ) and results in a 
positive brightness decline relation in this case. 

However, even an effective opacity declining with temperature may result in an anti-Amis- 
relation, because Amis depends both on the opacity effect and the excess energy available. For 
example, the pure detonation model DETl has an a ~ 0.73 and Amis of 1.37 mag. In contrast 



another pure detonation model DET2 has an a ~ 1.2 and Amis ~ 1-7 (jKhokhlov et al 



1993 



Hoflich &: Khokhlovlll996l ). In DETl, the high central density leads to significant electron capture 
in the central region producing a ^^Ni free core. In contrast, lower central densities in DET2 
produce ^^Ni in the center, and less ^^Ni in the outer region. In both models, the density structure 
is close to similar polytropes but the ^^Ni distribution is shifted outwards and inwards, respectively. 
The differences in a can be understood in terms of the deviation from Arnett's one-zone model 
with a flat ^^Ni distribution in mass. In DETl, this distribution is shifted outwards leading to an 
increased escape probability of 7-rays and greater expansion work which leads to reduced a. In 
DET2, the shift of the ^^Ni distribution inwards leads to a larger value of a. 

From the above discussion, it is obvious that a narrow Amis-relation requires similar abun- 
dance patterns and ^®Ni distributions for a given Mni. The theoretical relation Amis(MNi) depends 
on the explosion scenario. Self-similarity in the light curves requires a self-similar transformation 
between models of differing brightness. Within the delayed detonation scenario, M^i depends 
mostly on the pre-expansion during the deflagration phase, since if the density is too high neutron- 
rich iron-group elements will be produced. However, as the nickel mass drops out of the standard 
0.5-0.6 Mq range, the ^^Ni distribution shifts inwards wit h decreasing Mni. This produces a 
Amis(MNi)-relation which is in agreement with observations (IHofiich et al.ll2002l ). For core-normal 
SNe la, the temperature and, with it, the opacity remains high well after maximum light, leading 
to a "shallow" Amis(MNi). Below a certain M^i, the nickel is very centrally condensed (there is 
no nickel hole) and the opacity drops rapidly soon after maximum leading to a steep Amis(MNi)- 
relation. In the case of fast decliners ^^Ni is only produced during the deflagration phase. 
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3. Scenarios for SN 2001ay 

3.1. Models with an increased mass 

SN 2001ay was suggested to fall within the class of "Super-Chandra" mass SNe la. Ob- 
viously, SN 2001ay does not obey the standard Amis-relation. Within the framework of the 
empirical Amis-relation, SN 2001ay sho uld be brighter than ob served by roughly 1.0 magnitude 



which can be ruled out (see Fig. 10 of iKrisciunas et a 



tative super-Chaiidrasekhar SNe la, 2007if (IScalzo et al 



(Yamanaka et al 



200S 



Tanaka et al. 



201C 



201 ih. In fact, comp ared to the pu 



Silverman et al 



201C; Yuan et al 



2011 



comp i 
20ld) 



Taubenberger et al. 



and SN 2009dc 
esti- 



2011 



mates of Mboi at maximu m light differ by abou t 1 magnitude. Estimates for the brightness of 



SN 2001ay imply 0.5 M^i (jKrisciunas et al.ll201ll ). We need an increase of t he diffusion time scales 
by a factor of 2 from Amis of 1.25 for typical SNe la (jPhillips et al.lll999l ) to Amis = 0.68. The 
diffusion times scale as, t^is oc oc and, for homologous expansion, r oc t-\ If we assume Mch 
for normal SNe la, this would imply a progenitor mass of 2M0 which is super-Chandra and in fact 
typical of the mass suggested for mergers, but then one must explain why some super-Chandras 
are extra bright and SN 200 lay is not. 

A stronger argument against a high mass progenitor comes from the rise time to decline ratio. 
SN 2001ay rises quickly to maximum, whereas an increased time scale would unavoidably produce 
a slow rise to maximum light. 



3.2. The case for an Mch mass WD with a faster expansion rate 

Here, we suggest another scenario outside the "classical" regime for SNe la, but within 
the delayed-detonation scenario of Chandrasekhar mass white dwarf progenitors. To reproduce 
SN 2001ay, we require a fast rising light curve followed by a decline slower than the slope of the 
instantaneous radioactive decay. We can produce models with a < 1. This requires a model with 
^^Ni shifted to high velocities. Since eventually, the luminosity is given by the rate of instantaneous 
energy deposition, such a model will show a slow decline by construction. 



4. Results 

4.1. Explosion Models and Light Curves 

We describe how a change of parameters can transform the results from that given by normal 
LCs to that of SN 2001ay, within the same physical picture. 

Even in the absence of a well agreed upon model for SNe la most of the basic observational 
properties, that is, light curves and spectra, of core normal SNe la can be understood within the 
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framework of thermonuclear explosions of Chandrasekh ar mass white dwarfs (WD) and, in par- 
ticular, the delayed-detonation scenario (|Khokhlovlll989l ) which provides a natural explanation for 
the wide variety and range of the ^^Ni production. To first order, the amount of ^^Ni produced 
in the explosion depends on the pre-expansion of the WD during the deflagration phase, which 
in spherical models, can be parametrized conveniently by the density at which the transition be- 
tween deflagration to detonation occurs. For a recent review of the rela t ion b etween observational 
properties and the underlying physical model, see, for example, Hofiich (j2006l ). 



We base our analysis on a spherical explosion model within the framework of the pulsating 
delayed detonation scenario. Our goal is to construct a model with reduced a and short diffusion 
time scales. We seek to obtain our goals by increasing the expansion ratio and shifting of the ^^Ni 
distribution to higher velocity. 



4.1.1. Explosion Models 



The spherical explosion model has been constructed to allow fits of optical light curves and 
spectra of SN 2001ay. 

Within the DD scenario, the free model parameters are: 1) The chemical structure of the 
exploding WD, 2) Its central density, pc at the time of the explosion, 3) The description of the 
deflagration front, and 4) the layer at which the transition from deflagration to detonation occurs. 

As reference, we started from the delayed detonation model 5p0y25z22.25, which has been 
found t o be a good starting point for core-normal SNe la, with r e spect to both spect ra and light 
curves (|Dommguez et al.1 12OO1I : iHofiichI I2OO2I : lOuimbv et aD bood : iHoflich et all I2OI0I 1 . This base 
model originates from a star with a main sequence mass of 5 Mq and solar metallicity. Through 
accretion, the C+0 core of the star has grown close to the Chandrasekhar limit. At the time of 
the explosion of the WD, its central density is 2.0x10^ g cm~^ and its mass is close to I.STMq. 
The transition density ptr has been identifie d as the main factor which determines the ^^Ni pro- 



duction and, thus, t he brightness of SNe la (IHofiichI 1 19951 : iHofiich et al 



1995 



Hofiich et al. 



2002 



Iwamoto et al.lll999l ). The transition density ptr from deflagration to detonation is 25 x 10^ g cm 



For SN 2001ay, we tuned the parameters. The reference model has been modified as follows: 
We reduced the central density, pc to 1 x 10^ g cm~^, which decreased the potential energy, as 
always, and thus, increased the radius of the WD by ~ 30% to 2320 km, and we increased the C/0 
ratio in the former He-burning core (M < 0.56 Mq) to 4. Both modifications increase the explosion 
energy to 1.69 x 10^^ ergs (1.69 foe), leading to a more rapid expanding envelope, increasing the 
rate of geometrical dilution, and shifting the ^^Ni to higher velocity. These effects are responsible 
for both the fast rise and slower decline of SN 2001ay. 



In classical DD models, the deflagration to detonation tra r isition, DDT, occurs in a WD 
already unbound during the deflagration phase ()Khokhlovl Il989l : lYamaoka et al.Ml992l ). In con- 
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trast, in pulsating delayed models, PDD, the WD remains bound after the deflagration phase, and 
the DDT occurs during pulsati on (Khokhlov 19931 : iHofiich et ah 1995). Other authors have sug- 



gested variations on the PDD (llvanov a et al 



1974 : iBravo Garcia-Senj 12009 : iBaron et alJl2008l : 



Bravo &: Garci'a-Senal2009l : lBravo et al...20 09l. but here when we discuss the PDD model, we specif- 
ically refer to the model of iKhokhlovl (jl993). Reducing the burning rate during the deflagration 
phase moves the model from the DD to the PDD regime in which the WD is bound at the end 
of the deflagration phase. We artificially reduced the Atwood number, which we approximate as a 
constant independent of density and composition, (see Appendix) from 0.2 to (0.14, 0.12, 0.10) for 
the series PDD_lla-c which all undergo weak pulsation — producing LCs with rise times between 
14 and 16 days with a slow decline from maximum. The best fit to the observations is given by 
PDD_llb. Its details are described below. While some of the observed trends that we find with the 
PDD model may be obtained in a 3-D DD model including the effects of rotation, exploring such 
models is beyond the scope of this work. We mention results from PDD_lla/c below when needed. 

The density and abundance structure of PDD_llb is given in Figs. [TH2J The pulsation leaves 
a shell of ~ O.OGMq of unburned C/0 (with C/0 ~ 1, only the region of central helium burn- 
ing has enhanced C/0), and O.SIM© of ^^Ni. For pulsating delayed detonat ion models, mixing 



during the pulsation may occur. In contrast to previous pulsating models ([Hoflich et al. 



1995 



Hoflich Khokhlovlll996l ). we didn't assume full mixing during the pulsation, but instead mixing 
was limited to one scale height from the position of the burning front during the pulsation. The 
mixed region between burned matter forms the layer where the detonation is born. As a result, 
the density and chemical structure are similar to spherical DD models and consistent with late 
time line profiles and SN-remnants. We abandoned full mixing in order to preserve the inner, 
unmixed region as required by late time line profiles and the imaging of the supernova remnant 
S Andromedae (jllofiich et al.l I2OO4I: iMotohara etaP bood : IPesen et al.l I2OO7I : lOerardv et aD I2OO7I : 
Maeda et al.ll2010l ). Our models leave some unburned C/0. All modifications combined result in 
broader lines. We note that our model is highly parametrized: Firstly, the deflagration phase de- 
pends critically on the initial condition and as discussed above, the problem o f how to suppress 
stron g mixing is still unsolved, though it may be related to high magnetic fields (jPennev &: Hoflich 
201 ll ). Second, the amount of mixing will depend on Rayleigh- Taylor (RT) instabilities during the 
pulsation and possibly, rotation induced sheer instabilities. 



4-1.2. Light Curves 

Let us examine the formation of the light curves and the comparison with SN 2001ay. First, 
consider the energy input by 7-rays, positrons, and adiabatic cooling due to expansion. We contrast 
PDD_llb with our reference model. Both have a similar structure with respect to the chemical 
layering and, thus, similar opacities. Due to the higher explosion energy, the expansion rate of the 
inner layers is larger by about 25-30%. PDD_llb shows an increased energy loss due to expansion 
work. The opacities for 7's depend only on the mass column height, rjv/, and the electron/nucleon 
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ratio which is close to 0.5. tm scales as v x resulting in a 50% lower optical depth due to 
the higher velocity. The normalized energy deposition and escape fractions of 7-rays are shown in 
Fig. [3l The high expansion rate results in strong heating by 7-rays in the outer region and a higher 
escape probability for 7-rays. Between 10 and 25 days, PDD_llb has a an escape fraction larger 
than that of the reference model by about a factor of 2. 

In Fig. m we give a comparison between light curves of PDD_llb and the reference model. 
PDD_llb rises faster and is brighter at early times. The higher escape of 7-rays and the increased 
expansion work leads to a < 1 and therefore, a lower peak relative to the rate of instantaneous 
7-ray input. The second ingredient reducing the value of a is the opacity difference between optical 
photons and 7-rays. They are 0.1 — 0.2 and 1/35 cm^ g~^, respectively. 7-rays leak into the center 
where the energy is trapped without contributing to the light curve at maximum (Fig. [3]). 

We therefore have a deficit in luminosity (with respect to the instantaneous 7-ray deposition 
rate), and the optical light curves approaches "from below". The same opacity mechanism 
responsible for the regular brightness-decline relation results in an "anti-Phillips" -relation! More 
^^Ni and with it, larger opacities will lead to steeper declines. 



4.1.3. Comparison to SN 2001ay 

In Fig. O the theoretical LCs in B and V are compared with the observations. The agreement 
is reasonable, and they meet the brightness limit imposed by the early non-detection in R. The 
theoretical B and V have been corrected for the redshift z of the host galaxy. As discussed in the 
introduction, the distance modulus of the host galaxy i s 35.5 ±0.1 mag, and the galactic foreground 



extinction E{B — V) = 0.026 mag. iKrisciunas et al.l (|201lh use a reddening of the host galaxy of 



E{B — V) = 0.072 mag with an Ry of 3.1 and 2.4 for the Galaxy and host galaxy, respectively, 
giving Ay = 0.253 mag. Using our theoretical color, {B — y)max = 0.0 mag, and from an an 
optimized fit, we find the host galaxy reddening E{B — V) = 0.02 mag with a global Ry = 3.1 for 
a total extinction Ay = 0.144 mag. 



4.2. Spectral Analysis 

Based on the explosion models, 7-ray transport, and light curves, we analyzed the spectrum 
of SN 2001ay at maximum light using the NLTE-code PHOENIX (see Appendix). 

We use the density, abundance structure, and 7-ray and positron deposition given by the 
explosion models and light curve calculations. For the reddening and distance modulus, we use the 
values found from our LC fit. For the redshift, we use z = 0.030244. Converged models required 
256 optical depth points. 

For a consistent absolute My magnitude of —19.1 mag, the resulting synthetic spectrum is 
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compared to the spectrum observed at maximum light in Figure [H The continuum colors are 
well reproduced. The synthetic value of B — V = -0.06 mag, roughly consistent with mag of 
SN 2001ay. In o rder to avoid the complications of two dust populations that was assumed by 



Krisciunas et al.l (|201ll ). we have taken E(B-V) = 0.096 and Ry = 3.1, so our reddening is a bit 
"bluer" that was assumed in the light curve comparison, but it is well within the 0.1 magnitude 
error of the photometry. 

The spectrum is dominated by single ionized lines of S II, Si II, Ca II, Fe II, Co II as well 
as blends of doubly ionized species in the blue. The synthetic and observed spectra show good 
agreement. The Doppler shifts of lines from elements undergoing incomplete oxygen burning include 
Si II A6355, Si II A5972 and S II A5468 and A5654, and the Ca II H+K and the IR triplet are 
reasonably well reproduced within 10, 20, and 30 A, respectively, which corresponds to a velocity 

shift of 500 1000 km s~^ at the measured velocity of 14,400 km s~^. The strength of the 

absorption components agree well. A large number of weak lines form the quasi-continuum in 
the near IR. The feature due to blends of Fe II/Co II/Ni II which could possibly be misidentified 
as due to C II 6580 is rea sonably strong, but the feature does not agree with the absorption 



notch identified as C I I bv iKrisciunas et al.l (j201ll ). However, the evidence for C II claimed in 



Krisciunas et al.l (j201ll ) is stronger than warranted based on the analysis that were done for that 
paper (R. Thomas, personal communication). Some disagreement is evident. The S II W at about 
5000 A, blended with other lines, and the Si II A5970 line are too weak. In the red, the Ca IR 
triplet is too weak, as is the O I AA7773 line, which is however, severely contaminated by Telluric 
absorption. The lack of O I is an indication that the outer part of our model is not correct, either 
in density structure or due to temperature effects. 

To probe importance of consistency and the sensitivity of the spectra, we calculated a series 
of maximum light spectra which are under- and over-luminous, with My = (-18.92,-19.23,-19.35,- 
19.45) mag. 

At +0.1 mag (Fig. [7l the B — V color becomes 0.03 mag. The Si II A6355 line is well-fit and 
the Ca II IR triplet is relatively well reproduced. In the blue, Co II lines dominate and a Co II 
line at 7400 A appears. The feature at 6500 A is weaker. The Doppler shift of the Si II and S II 
remain in good agreement because they are formed within the region of incomplete oxygen burning 
which produced nearly constant abundances due to burning in quasi-statistical equilibrium (QSE) 
conditions (Fig. [2]) . 

At -0.15 mag (Fig. B — V equals —0.11 mag, but is clearly a little bit too blue as can be 
seen in the red part of the spectrum However, the Si II A5970 feature is too weak and S II lines 
are still significantly too weak. The Ca II IR triplet is much too weak. In the blue, Fe II lines 
are stronger and the fit is somewhat better. This could imply that the model should have had an 
initial metallicity somewhat higher, which would create more Fe II at the expense of Co II (from 
radioactive nickel). 

At -0.25 mag (Fig. [9]), the spectrum becomes too blue, B — V = -0.08 mag, and the UV flux 
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becomes too large. The Si II and S II lines are still significantly too weak and the IR flux is too 
low. The Ca H+K line is now too weak. 

At -0.35 mag (Fig. [TOj) . the spectrum is much too blue, B — V = -0.07 mag and the Ca II 
features are completely absent. The Si II and S II features are now much too weak. 

Figure [TT] shows the model fit to the observed HST +Keck spectrum of April 29. The Si II line 
is very nicely fit, the S II "W" is well fit, although the blue line is a bit stronger in the model than 
in the observed spectrum. The rest of the blue is reasonably well fit, although the flux is about 
a factor of 2 too high in the UV, which may be an indication of the need for a higher metallicity, 
which would produce more line blanketing and reduce the UV flux. Both the Ca H+K and IR 
triplet are reasonably well flt, and the Co II + Fe II + Mn II feature at about 7350 A is a bit 
too strong. The O I AA7773 line is much too weak. Figure [12] shows the model flt for a somewhat 
brighter model My = —19.22 and B — V = 0.22, still a bit too red, but the line features flt similarly 
well to the redder case. Figure [T3] shows the model flt for a somewhat brighter model My = —19.40 
and B — = 0.10, a bit too bright, but about the right color. Again, the too high UV and the 
weak O I line are indications that the outer structure of our model is not quite correct, which could 
be due to primordial metallicity or our assumptions about mixing in the FDD. 

Obviously, spectral synthesis is a very sensitive tool to probe the structure and luminosity on 
an level of AMy ~ 0.1 mag. Variations in excess of AMy ~ 0.2 lead to spectral flts which would 
lead to a rejection of even a valid, underlying explosion model. We note that our best flts give a 
"correct" velocity shift of e.g. the Si II feature because, in PDD_llb, it is formed in the region of 
QSE for the Si/S group. Consistent treatment of explosion, light curves, and spectral formation is 
important. However, spectra lines are sensitive to both temperature effects and line blending and 
even a strong line like Si II A6355 is not unblended. One can see a small blueward shift between 
Figs. [6] and [TTJ Figure [U] shows the Si II A6355 feature in velocity space, where two models at 
the same epoch (Apr 29) are plotted along with the data. The cooler model has an absorption 
minimum that is bluer than that of the hotter model and both are a bit bluer than the observed 
spectrum, although the noise in the spectrum makes a determination of the minimum a bit difflcult. 

In turn, light curves are mostly sensitive to the inner region and the large number of weak 
lines. Both in the LC and spectral calculation, B — V '\s consistent, namely and -0.06 mag. Only 
spectra show that the Si II line is too weak. 

As discussed above, spherical models inherently suppress mixing by instabilities due to the 
deflagration phase at the inner layers, and the interaction between shells and the outer layers 
during the acceleration phase. 

Spectral analysis requires underlying models which are consistent including the luminosity. 
Then, however, they allow the study of secondary effects of individual spectral features. 
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5. Discussion and Conclusion 

The peculiar Type la SN 2001ay is an important milestone for our understanding of the 
explosion, light curves, and spectra. 

Even with its extremely unusual light curve shape and spectral features, we showed that it 
still can be accommodated within the framework of the explosion of a Chandrasekhar mass white 
dwarf. 

Although it does not follow the Amis-relation, the light curve can be reproduced within the 
physics of normal SNe la, and in the framework of pulsating delayed detonation models. In our 
models, the key difference is a high, 80% carbon mass fraction, rather than the 15-20% carbon 
mass fraction usual for stellar central He-burning. The size of the carbon rich core is 0.56 Mq. The 
excess of carbon coupled with a lower central density of the initial WD results in an increase of the 
final kinetic energy by about 40% and in turn, shifts the distribution of ^^Ni to larger velocities. 
These modifications are responsible for both the LC characteristics and the broad spectral features. 

The peculiar light curve shape is a consequence of the shift of the ^^Ni distribution in velocity 
space and higher rate of central expansion. Transport effects of 7-rays become more important and 
the escape probability is increased. This leads to a fast, early rise and a value of a < 1, where 
a is the ratio between the luminosity and the instantaneous energy generation rate by radioactive 
decays at maximum. In normal SNe la, the Amis-relation can be understood as a consequence of 
the temperature dependence of the opacity in combination with a > 1. As discussed in § [21 more 
^^Ni causes higher luminosities and temperatures. Higher temperatures lead to larger opacities 
and therefore, a slower drop of the luminosity from the typical value of ~ 1.2x the instantaneous 
energy input rate at maximum, to the instantaneous energy input rate at later times. Accordingly, 
for SN 2001ay with a ~ 0.55, the light curve approaches the instantaneous energy input rate from 
below which explains the unusually slow decline. Within our models and as a corollary, we expect 
that there exists a sub-class of supernovae which obeys an anti-Phillips relation. As discussed 
below, this sub-class should be rare. 

The unusually broad spectral features can be understood by the overall shift of the overall 
element pattern to higher velocities. However, this pattern remains similar to that found in normal 
SNe la (Fig. [2]). We demonstrated the power and sensitivity of spectral analysis. High sensitivity 
demands a consistent treatment with the explosion and light curve models. Otherwise, even valid 
models may be rejected because of poor fits. However, since a good overall fit is achieved with 
a consistent treatment, individual features are a powerful tool to study details of the explosion 
physics. 

We have shown that limiting the mixing during the pulsation to small scales produces a very 
similar abundance pattern to the one produced in standard DD-models which have been shown to 
reproduce the observables for the majority of SNe la. PDD may be a promising mechanism for the 
DDT. 
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While our model fits the basic observed features of SN 2001ay, we do not mean to imply that 
other models are not possible. In particular a DD model including 3-D effects as discussed above 
with a high central C/0 ratio may also be a viable model and there may exist other models that 
fit the observed trends. Pursuing other explanations is beyond the scope of the present work. 

Our model agrees reasonably well with the observations (Fig. [5]) but the high carbon abundance 
in the progenitor poses a challenge. Th e central region of th e progenitor originates from central He 
burning in stars with less than 7-8 Mq (iBecker &: IbenlllOSOl ) . During the early stages of central He 
burning, high carbon abundances are produced by ^He(2Q!, 7)^'^C burning. With time, the helium 
abundance is reduced in the core. Th en, ^'^ Cja^iY^O depletes to a typical value of 10-25% 
(jUmeda et al.lll999l : iDomfnguez et ahlboOlh . The final amouri t of depends on the mass 

of the pr ogenitor, and the ^"^Cia.^Y^O rate ('Buchmann 1997 ), and the amount chemical mixing 



2001 



assumed (Castellani et al. 



1985 



Caputo et al.i.l989 : .Renzini fc Fusi Pecci 



1988 



Domfnguez et al 



Increased chemical overshooting (or semiconvection) prolongs the phase of burning un- 
der helium depleted conditions, leading to a lower value of C/0. Although chemical overshooting 
prescriptions vary widely between various groups studying stellar evolution, the final outcome is 
a C-poor mixture. To reac h high central carbon abundances, burning un der He-depleted condi- 
tions must be avoided (see IStraniero et al.l l2003l . and references therein) . IStraniero et al.l (j2003l ) 
found that they could increase the central carbon abundance somewhat by increasing mechani- 
cal overshooting. As a possible solution, we propose a common envelope phase with very strong 
mixing induced by a compact secondary such as a brown dwarf or planet. Common envelope 
evolution is generally assumed to be responsible for the forma tion of close binaries leading to cata- 
clysin ic variables, X-ray binaries, and supersoft-X-ray sources (iLivio &: Sokerlll988l : iRicker fc: Taam 
20081). Moreover, close biriary systems with plan et or brown dwarfs have recently been detected 



(jNeuhauser k. Posseltl 120071 : iHessman et al.l l201ll ) . Unfortunately, there are no systematic studies 



which allow one to estimate the amount of mixing. Detailed, numeric al studies are few, and ana- 
lytic models are insufficient to quantify the amount of chemical mixing (jMever &: Mever-Hofmeister 
19791 ). One observ ational clue may be SN 19 87A, which is believed to be the result of a common 



licity (iBrunish fc Truran 



1982 



Chieffi et al 



envelo pe evolution (IPodsiadlowski et al.lll99C ). Blu e progenitors may also be a result of low metal- 



2003), and very few SN1987A-like events have been 



found (|Pastorello et al.ll2012l : iTaddia et al.ll2012l ). iPastorello et al.l (|2012l ) estimate that SN 1987A- 
likes represent ~ 1.5 — 4% of SNe II; however, this is likely an extreme upper limit to our scenario, 
since the preponderance of these objects may be just from higher mass compact stars, rather than 
from systems in binaries that have undergone common envelope evolution. Similarly, one can use 
observations of white dwarfs to provide anothe r estimate. The number of white dwarfs with sub- 
stellar companions is < 0.5% ()Farihi et al.ll2005l ), and the number o f white dwarfs with debris disks 
from tidally disrupting minor planets is 1 — 3% (jFarihi et al.ll2009l ) Therefore, our scenario should 
be quite rare, representing 0.05 — 0.5% of all SNe la. 



Finally, we want to mention the long list of limitations of our studies: We have suggested 
the existence of a rare sub-class of SNe la which should obey an anti-Phillips relation. It is up to 
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future, systematic surveys, such as LSST to find a sufficiently large statistical sample. Note, that 
this sub-class may be hard to separate at the bright end of SNe la because the decline rates for 
both normal and "SN 2001ay-like" supernovae are low and similar. Chemical mixing during the 
common envelope phase needs to be studied in detail. Our model for SN 2001ay requires a shift of 
the ^^Ni distribution and a central hole in ^^Ni. Although verified for a number of normal SNe la, 
late time spectra for a SN 2001ay-like are required to confirm our assumptions of little mixing. 
Detailed timc-scrics of early time spectra may help to probe whether PDDs provide a common 
mechanism for the transition from deflagration to detonation. Detailed 3D models of the pulsation 
phase are required to test and quantify possible mixing. 

The work presented in this paper has been carried out within the NSF project "Collabora- 
tive research: Three-Dimensional Simulations of Type la Supernovae: Constraining Models with 
Observations" whose goal is is to test and constrain the physics of supernovae by observations 
and improve SNe la as tools for high precision cosmology. The project involves The University 
of Chicago (AST-0709181), the University of Oklahoma (AST-0707704), Florida State University 
(AST-0708855), Texas A&M (AST-0708873), The University of Chile in Santiago, and the Las 
Campanas Observatory, Chile. This research was also supported, in part, by the NSF grant AST- 
0703902 to PAH. The work of EB was also supported in part by SFB 676, GRK 1354 from the DFG, 
and US DOE Grant DE-FG02-07ER41517. ID has been supported in part by the Spanish Min- 
istry of Science and Innovation project AYA2008-04211-C02-02 (ID). Support for Program number 
HST-GO-12298.05-A was provided by NASA through a grant from the Space Telescope Science 
Institute, which is operated by the Association of Universities for Research in Astronomy, Incor- 
porated, under NASA contract NAS5-26555. This research used resources of the National Energy 
Research Scientific Computing Center (NERSC), which is supported by the Office of Science of the 
U.S. Department of Energy under Contract No. DE-AC02-05CH11231; and the Hochstleistungs 
Rechenzentrum Nord (HLRN) . We thank both these institutions for a generous allocation of com- 
puter time. 
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A. Brief Description of numerical Methods 
A.l. Explosion 

We have calculated explosion models and light curves using the one- dimensional radiatio n- 
hydr o code HYDRA using computational modules for spherical geometry ( Hofiich 1995 . 2003al lbl. 
20091 ). We solve t he hydrodynamical equatio ns explicitly by the piecewise parabolic method on 
910 depth points (IColella Woodward! 1 19841 ). Because a simple a network is insufficient to de- 
scribe in sufficient detail the chemical boundary. Nuc lear burning is taken in t o account using an 
extended network of 218 isotopes from n, p to ^'^Kr ( Thielemann et al.lll996l : iHix fc Thielemann 
19961 : iHoflich et al.l Il998l . and references therein). The propagation of the nuclear burning front 



is given by the velocity of sound behin d the burning front in t he ca se of a detonation wave. We 
use the parametrization as described in iDommguez &: HofiichI (j200d ). For a deflagration front at 
distance rbum from the center, we assume that the burning velocity is given by Wburn = max(ft, f;), 
where vi and vt are the laminar and turbulent velocities with 



VargLf, 



[1] 



with 



and 



OT = {a — l)/{a + 1) 



a 



urn J ■ 



Here ot is the Atwood number, Lj is the characteristic length scale, and /o"*" and p~ are the densities 
in front of and behind the burning front, respectively. The quantities a and Lf are directly taken 
from the hydrodynamical model at the location of the burning front and we take Lf = rburn(i)- 
The transition density is treated as a free parameter. Table [T] gives the yields of the stable elements 
as well as the total amount of ^^Ni produced in the explosion. 



This preprint was prepared with the AAS IATJhjX macros v5.2. 
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Table 1. Element Yields 



Element 


Yield (M0) 


He 


0.01 


C 


0.021 





0.043 


Ne 


3 X 10~^ 


Na 


7.8 X 10"^ 


Mg 


1 X 10"^ 


Si 


0.182 


P 


9 X 10"^ 


S 


0.104 


CI 


4.2 X 10~^ 


Ar 


0.022 


K 


4.8 X 10-5 


Ca 


0.022 


Sc 


3.1 X 10"'^ 


Ti 


1.5 X 10-3 


V 


0.022 


Cr 


0.106 


Mn 


0.024 


Fe 


0.661 


Co 


1.1 X 10-3 


Ni 


0.146 


56^1 


0.515 



Note. " The yield 
of the stable elements at 
time infinity along with 
the total mass of ^^Ni 
produced in the explo- 
sion. 
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A. 2. Light Curves 



From these explosion models, the subsequent expansion, bolometric and broad band light 
curves (LC) are calculated using monochromatic radiation transport via the Eddington Tensor 
method. Both the Eddington tensor and the 7-ray transport are calculated via Monte Carlo as 
in the 3D case. We include scattering, free-free, and bound-free opacities, and include the line 
transitions in the Sobolev approximation. For several elements, including C, O, Mg, Si, S, Ca, Fe, 
Co, and Ni we solve the statistical equations for the three main ionization stages using detailed 
atomic models with 10-50 super-levels with ~ 12,000 transitions using the databases of iKurucz 



(|2002l ) and ICunto &: Mendozal (1 19921 ) . The levels closest to the continuum and of other elements 
are treated in local thermodynamical equilibrium with about 10^ line transitions. 



A. 3. Spectra 



Th e calculations were perforrned using the multi-purpose stellar atmospheres program PH DENIX/ ID ver- 
sion 16 (jHauschildt BaronlE999l : baron HauschildtlEggsl : iHauschildt et al.lE997al lbl .ll99d). Ver- 



sion 1 6 inco rporates many changes over previous versions used for supernova modeling (jBaron et al 



20071 . 12OO6I ) including many more species in the equation of state (83 versus 40), twice as many 
atomic lines, and many more species treated in full NLTE and an improved equation of state. 
PHOENIX/ID solves the radiative transfer equation along characteristic rays in spherical symmetry 
including all special relativistic effects. The non-LTE (NLTE) rate equations for many ionization 
states are solved including the effects of ionization due to non-thermal electrons from the 7-rays 
produced by the radiative decay of ^^Ni, which is produced in the supernova explosion. The atoms 
and ions calculated in NLTE are: He HI, C I-III, O I-III, Ne I, Na I-II, Mg I-III, Si I-III, S I IU, 
Ca II, Ti II, Cr I-III, Mn I-III, Fe I-III, Co I-III, and Ni I-III. These are all the elements whose 
features make important contributions to the observed spectral features in SNe la. 7-ray deposition 
was that calculated by the light curve. 

Each model atom includes primary NLTE transitions, which are used to calculate the level 
populations and opacity, and weaker secondary LTE transitions which are included in the opacity 
and implicitly affect the ra te equations via their effect on the solution to the transport equation 
(jHauschildt &: Baronlll999l ). In addition to the NLTE transitions, all other LTE line opacities for 
atomic species not treated in NLTE are treat ed with the equival ent two-level atom source function, 
using a thermalization parameter, a = 0.10 (jBaron et al.lll996l ). The atmospheres are iterated to 
energy balance in the co-moving frame; while we neglect the explicit effects of time dependence in 
the radiation transport equation, we do implicitly include these effects, via explicitly including p 
work and the rate of gamma-ray deposition in the generalized equation of radiative equilibrium and 
in the rate equations for the NLTE populations. 

The outer boundary condition is the total bolometric luminosity in the observer's frame. The 
inner boundary condition is that the flux at the innermost zone {v = 700 km s-i) is given by the 
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diffusion equation. Converged models required 256 optical depth points to correctly obtain the 
Si II A6355 profile. 

1995al lbl. 



P HOENIX and HYDRA have been well tested and compared on SNe la (iNugent et al 



19971; iLentz et al.l 120011 : Baron et al.ll200fil : lHofiichll200,4 iHofiich et al.lll998l : iHoflichI I2OO2I I and, in 



particular compared with observed light curves and spectra of SN 1994D and SN 2001ay. 
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Fig. 2. — Abundances of stable isotopes as a function of the expansion velocity for model PDD_llb. 
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Fig. 3. — We show the normahzed energy deposition by 7-ray and positrons in PDD_llb (left), 
and the escape fraction for PDD_llb and the reference model (right). 
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Fig. 4. — We give the instantaneous deposition rate due to radioactive decay (dashed) and bolo- 
metric luminosity (sohd) for PDD_llb, a pulsating model with a carbon rich core, and a "classical" 
DD model. The light curve of the DD model has been uniformly shifted by -0.037 dex and the light 
curve and 7-ray deposition of the PDD_llb model have been uniformily shifted by +0.05 dex to 
better illustrate the variation of the gradients. 
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Fig. 5. — B and V LCs of SN 2001ay (jKrisciunas et alJl201ll ) in comparison with theory. The 
comparison between SN 2001ay and PDD_llb in B (upper) and V (lower) as a function of time 
since maximum hght in the F-band. The model light curves have been cor rected for reddeii i ng ari d 
redshift. We assume a distance modulus m-M of 35.55 mag, and following iKrisciunas et al.l (|201ll ). 
we take Ay = 0.253 and Ab = 0.35. The data includes the photometric errors. 
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Fig. 6. — Synthetic spectrum at day 16 of PDD_llb in comparison with maximum spectrum of 
SN 2001ay. The density and abundance structure and the energy input by radioactive decay are 
taken from the explosion and hght curve calculation (Figs. [2l [3l and[5|). Roughly consistent with 
the light curve, My was taken to be -19.07 mag, we have assumed standard reddening Ry = 3.1 
and E(B-V) = 0.096. 
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Fig. 7.— 



Same as Fig. [6]but dimmer by 0.1 mag, i.e. My of -18.92 mag. 
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Fig. 9.— 



Same as Fig. [6] but brighter by 0.25 mag, i.e. My of -19.35 mag. 
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Fig. 11. — Synthetic spectrum (red) for 23 days after explosion compared with the observed spec- 
trum obtained on Apr 29, 2001, from both HST and Keck. My is -18.97 mag. and B - V = 0.43. 



- 33 - 




- 34 - 




le-15 
1.6F 

1.4 




Velocity (km s ^ ) 



Fig. 14. — The Si II A6355 feature for 23 days after explosion compared with the observed spectrum 
obtained on Apr 29, 2001 (black) with the model shown in Fig. [13] (green) and a somewhat hotter 
model (blue). The blueward shift in the absorption is clearly a temperature effect and is also due 
to blending. 



